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ABSTRACT 

Carbon-rich grains with isotopic anomalies compared to the Sun are found in primi¬ 
tive meteorites. They were made by stars, and carry the original stellar nucleosynthesis 
signature. Silicon carbide grains of Type X and C, and low-density graphites condensed 
in the ejecta of core-collapse supernovae. We present a new set of models for the explo¬ 
sive He shell and compare them with the grains showing and ratios 

lower than solar. In the stellar progenitor H was ingested into the He shell and not 
fully destroyed before the explosion. Different explosion energies and H concentrations 
are considered. If the SN shock hits the He-shell region with some H still present, the 
models can reproduce the C and N isotopic signatures in C-rich grains. Hot-CNO cycle 
isotopic signatures are obtained, including a large production of and The 
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short-lived radionuclides and ^®A1 are increased by orders of magnitude. The pro¬ 
duction of radiogenic ^^Ne from the decay of ^^Na in the He shell might solve the puzzle 
of the Ne-E(L) component in low-density graphite grains. This scenario is attractive for 
the SiC grains of type AB with ratios lower than solar, and provides an alter¬ 

native solution for SiC grains originally classified as nova grains. Finally, this process 
may contribute to the production of and in the Galaxy, helping to produce the 
ratio in the solar system. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: abundances 
— stars: evolution — stars: interiors — supernovae: general 


Introduction 


Primitive meteorites are carriers of several types of dust of presolar origin, coming from different 
stellar sources. These sources are identified through measurements of the isotopic abundances in 


single grains (Zinner 2014). Presolar C-rich grains that condensed in Core-Collapse Supernovae 


(CCSNe), are low-density (LD) graphite grains (carrier of the Ne-E(L) component, Amari et al. 


1990), nano-diamonds (carrier of the Xe-HL component, Lewis et al. 1987), and silicon carbides 


(SiC) of Type X (about 1% of all presolar SiC grains, e.g., Besmehn Sz Hoppe 2003) and of Type 
C (0.1-0.2% of all presolar SiC grains, e.g., Pignatari et al. 2013b). Evidence for initial “^^Ti (a 


prominent signature of SN grains) led to the conclusion that at least one of the few SiC grains 


originally classihed as nova grains is made in CCSNe (. 

Mittler & Hoppe 

2005 

). The SiC grains 

of Type AB (about 4-5% of all presolar SiC grains. 

Amari et al. 

2001 

possibly have multiple 

astrophysical sources, including CCSNe (e.g., 

Zinner 

201^ 

1). 




15 


N-excesses are a signature of SiC grains from CCSNe. They are also found together with 


18 


O excesses as hotspots, probably linked to TiC subgrains, in LD graphite grains (Groopman 


et al. 2012). The combination of both and ^^C/^^C ratios lower than solar observed i 


m 


a fraction of C-rich grains is a challenge to theoretical CCSNe models, assuming mixing between 


different CCSNe layers and maintaining C/0>1 (e.g., Travaglio et al. 1999). ^®N production i 


m 


the explosive He-burning shell does not provide a satisfactory solution (e.g., Bojazi &: Meyer||2014 ). 
Presolar SiC grains originally classified as nova grains have ^^N-excesses higher than those of SN 
grains, in accordance with predictions from nova models. In contrast, a satisfactory explanation for 
the origin of ^®N-excesses in some SiC-AB grains has not been found yet (Huss et al.||l997 Zinner 


2014). 


The requirement of a C-rich environment to condense C-rich grains poses a serious problem 


for finding scenarios that explain the observed isotopic abundances in SN grains. Yoshida (2007) 


reproduced the C and N isotopic ratios of a large number of SiC-X grains by mixing matter from 
different SN layers. Although these mixtures have C/0>1 in many cases, the C>0 constraint and 
the required heterogeneous and deep mixing in SN ejecta is problematic. To overcome this issue, 
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(2010) showed that the isotope ratios of SiC SN grains are not consistent with predictions for 0-rich 
matter. 


In this work we show that the entrainment of H-rich material into the He shell before the 
SN explosion allows the coproduction of and ^®A1, accounting for isotopic signatures of 

presolar SiC grains with low ratios. This also provides a new production scenario for SiC 

grains classified as nova and AB grains. 

The paper is organized as follows. In ^ we describe stellar models and nucleosynthesis cal¬ 
culations, in ^ theoretical results are compared with isotope data for different types of presolar 
grains. In ^ we discuss the implications of these calculations for the nitrogen inventory in the 
Galaxy. In ^results are summarized. 


2. Stellar model calculations and nucleosynthesis 


A set of CCSN nucleosynthesis models is calculated from a progenitor 25Mq, Z = 0.02 star 
( ]Pignatari et al.|2013c P13 hereafter). The simulation of the CCSN explosion includes the fallback 
prescription by Fryer et al. (|2012 ) (model 25d). The initial shock velocity beyond fallback is 


Vg = 2 X 10®cm s“^, which corresponds to a total explosion energy of Egxp = 4 — 7 x lO^^erg 
After the end of the central C burning, the convective He shell becomes unstable and eventually H 
from above the shell is ingested into the He-rich region. The onset of core O burning completely 
deactivates the convective shell until the SN explosion, leaving He-rich shell material with about 
1.2% of H. The occurrence of mergers between different stellar zones in massive stars is known 


from the literature (e.g., Woosley & Weaver 1995, for H mixed into the He shell), but its effect on 


pre- and post-explosive nucleosynthesis has never been investigated before. The ingestion of H into 
hotter layers needs to be simulated with multi-dimensional hydrodynamics because one-dimensional 
hydrostatic models lose their predictive power ( Herwig et al.|[Ml4 ). In the absence of realistic 3D- 
hydrodynamics simulations, a qualitative analysis of the effect of residual H during explosive He 
burning is required. Together with model 25d from P13 we consider five new models: 25d-H5, 
25d-H10, 25d-H20, 25d-H50 and 25d-H500 (models set d). These models share the structure and 
CCSN explosion of 25d, but the amount of H in the He shell is reduced by a factor of 5, 10, 20 
50 and 500 respectively compared to model 25d. In P13, the CCSN model with the most extreme 
conditions is the 15Mq, Z = 0.02 model, with rapid SN explosion. The bottom of the C-rich 
C/Si zone reaches a peak temperature 3.3 times higher than that in the 25d model (2.3xlO®K 
compared to 0.7xl0®K). The peak density is 100 times larger. We produced a new set of 25Mq 
models, artificially increasing temperature and density to mimic the He-shell conditions of the 
15Mq model, and exploring the same range of pre-supernova H concentrations in He-rich material: 
25T-H, 25T-H5, 25T-H10, 25T-H20, 25T-H50 and 25T-H500 (models set T, see Table[^. 
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The post-SN abundances of H, He, C and N isotopes, and of the short-lived 

species ^^Na, ^®A1, “^^Ti and ^'^Fe are reported in Fig.[^ for the models 25T-H and 25T-H50. 
Model 25T-H has the largest H concentration and the highest explosion temperature and density. 
At the bottom of the He shell at ~6.8Mq, ^®Si is made by a— capture on forming a C/Si 
zone (Pignatari et ah 2013a). Between 6.83 and 7 .O 4 M 0 the abundance of drops, forming an 


0-rich zone. We call this the 0/nova zone. It has C/0<1, and the abundances carry an explosive 
H-burning signature. Qualitatively there are features similar to those of nova nucleosynthesis 


(e.g., Jose & Hernanz 2007) except in the deepest region of the 0/nova zone, which shows unique 
nucleosynthesis features. In particular, ^^C(p, 7 )^^N depletes while feeds ^^O and ^®0 via 
the ^^N(p, 7)^^0 and the ^^N(a,p )^®0 reactions. 

Similar nucleosynthesis of explosive H burning is obtained in the 25d models but within a much 
smaller region (<O. 1 M 0 ), and most of the He-shell material is not affected by the explosion. These 
layers are mostly carrying the pre-SN signature of H burning under He-shell conditions. The 
and abundances are 0.5% and 0.3% respectively. A mass fraction of a few 10“® and 10“^ is 
obtained for ^®A 1 and ^^Na. 


In the model 25T-H protons suppress the SN neutron burst in the He shell (e.g., Meyer et al 


2000). Protons power the ^^Ne(p, 7 )^^Na reaction, in competition with the ^^Ne(a,n)^^Mg reaction. 
The presence of pre-explosive provides the ^^C(a,n)^®0 alternative neutron source, mitigated 
by the presence of the neutron poison ^^N(n,p)^‘^C, and by the competition of the ^^C(p, 7 )^^C 
reaction. The H-burning products ^^Na and ^®A1 are boosted (and X(^'^Ti)<10“® is made from 
proton captures on Ca isotopes), while neutron-burst products such as ®°Fe are reduced (Fig.[^ 
bottom panel). This effect depends on the explosion energy and on the H abundance. In model 
25T-H the ^®A1 peak is three orders of magnitude larger than in model 25T-H50. The ^^Na peak 
is about hve orders of magnitude larger. On the other hand, ®*^Fe is efficiently produced only in 
model 25T-H50. These results might have implications for the observation of the 7 —emission from 
the decay of ^®A1 and ®^Fe in the Galaxy and for their relative abundances ( Diehl|[2006 ). 

In our simulations ^^Na is made in the 0/nova zone by proton captures, with ^*^Ne being 
made by a-capture on ^®0 as the main seed. In model 25T-H the abundance peak for ^^Na is 


almost 10%.This scenario may explain the Ne-E(L) component in LD graphite grains (Amari et al. 
1990). This component consists mostly of ^^Ne, suggesting that it is of radiogenic origin from ^^Na 


originally condensed into the grains. In baseline massive star models, ^^Na is made by C burning in 
the deeper 0-rich 0/Ne zone, however it was unclear how to mix Na into the C-rich zones without 
admixing large quantities of O. For the first time, we identihed a potential pathway to produce 
^^Na in the He-shell layers. 
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3. Comparison with presolar grains 


In this section, the nucleosynthesis models described in ^are compared with measurements of 
different types of C-rich grains with ratios lower than solarfrom the Washington University 


Presolar Grains Database (Hynes &; Gyngard 2009 and original references in the database). We 


focus on the abundances in the G-rich regions in the He shell (the C/Si zone and the C-rich He/C 
zone above) and in the 0/nova zone, which lies in between. 

In Fig.[^ C and N isotopic ratios of presolar SiC grains are shown together with the predictions 
from our models in the He shell (He/C zone). For ^^C/^^C lower than solar, the theoretical curves 
for the models 25T-H5, 25T-H10, 25T-H20, 25d-H5, 25d-HI0, and 25d-H20 reproduce the ratios 
observed in nova grains, and are compatible with X grains with low N isotopic ratios assuming 
some degree of dilution with material of close-to-normal composition. The 25T-H model produces 
the lowest ^^C/^^C ratios, due to large contributions from radiogenic ^^C from (Fig.[^. The 
25T-H5, 25T-H10, 25T-H20, 25d-H5, 25d-H10, and 25d-H20 models show similar ^^C/^^C ratios. 
The lower limit of the ratio depends on both explosion energy and H concentration. 

requires high temperatures to be made as a radiogenic product of The 25d-H10, 25d-H20, 
25d-H50 and 25d-H500 models directly reproduce for the first time the isotopic compositions of 
SiC-AB grains with N isotopic ratios lower than solar. The 25d and 25d-H5 models show ratios 
lower than observed, and dilution with normal material is needed. Most of graphites show the 
well-known signature of N isotopic equilibration or contamination (Zinner 2014[ ). Our predictions 
for C and N ratios might also provide a solution for mainstream SiC grains with ^^N/^®N<100, not 
reproduced from low-mass Asymptotic Giant Branch star models ( Palmerini et al.||2011 ). Whether 
these grains could have condensed from SN ejecta needs to be explored. 

In Fig.j^ upper-left panel, ^^C/^^C and ^®A1/^^A1 ratios of C-rich grains are compared with 
the ratios from the 25T models. For our models the curves range from high ^^C/^^C ratios in the 
C-rich C/Si zone to low ^^C/^^C ratios in the C-rich He/C zones, with the 0/nova zone having 
12 c/13c ratios in between. Therefore, C-rich mixtures from the C/Si, 0/nova, and He/C zones of 
the models 25T-H, 25T-H5, and 25T-H10 can potentially reproduce for the fist time the highest 
^®A1/^’^A1 ratios observed in SiC-X and LD-graphite grains, consistently with the observed 
ratio. In the upper-right panel we compare model predictions with LD graphites for ^®0/^®0 and 
12 c/13c_ also report the results for the ISM© model (Fig. 3, Pignatari et al. 2013a), which 
does not include H-ingestion and has CCSN conditions similar to those of the 25T models (® . As 


discussed in Pignatari et al. (|2013a|), in the most ^^C-rich layers there is negligible amount of oxygen 


left and minor mixing with normal material can yield a close-to-normal O isotopic composition. 
The effect of the H ingestion is to destroy and increase ^^C, moving the theoretical curves for 
the He/C zone from the lower-right quadrant to the upper-right. Unfortunately, a large fraction 
of graphites show O isotopic equilibration or contamination, in particular for graphites with the 


lowest ^^C/^^C (e.g., Zinner 2014), making a meaningful comparison between grain data and models 
difficult. 
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The 25d models seem to be more compatible with SiC-AB grains (Fig.j^. In Fig.j^ lower 
panels, we compare isotopic ratios of AB and nova SiC grains with theoretical curves for 25d 
models. AB grains with lower than solar are on average more ^®Al-rich. A nova grain 

and the AB grains with the highest ^®A1 enrichment are compatible with our models, while most of 
them require some dilution with material of close-to-normal composition. This is consistent with 
an observed trend of ^®A1/^^A1 inversely proportional to the A1 abundance. Theoretical predictions 
for ^®Si/^®Si and ^^Si/^®Si (not shown in the figure) and C isotopic ratios seem to be consistent 
with AB and nova grains. 


Previously applied SN mixing schemes considered matter from the He/N zone to account for 
low together with high ^®A1/^^A1 in SN grains (e.g., Zinner 2014). However, in contrast 

to our 25T and 25d models, those calculations do not provide in sufficient quantities. The 
possibility of having large concentrations of and ^®A1 in the He/C zone makes it easier to 

fit the observed ratios within the constraint of a C-rich mixture in mixing models. 


4. Production of nitrogen in the Galaxy 


The peculiar isotopic abundances discussed in the previous sections are relevant for galactic 
chemical evolution (GCE). While presolar SN grains found in meteorites only condensed from 


CCSNe that exploded at most several hundred million years before the formation of the Sun (Jones 


et al. 1994), H-ingestion events in the He shell are expected to happen more frequently in low- 


metallicity stars. About 15% of SiC-X grains and 53% of LD graphite grains show lower than solar 
C isotopic ratios (see Presolar grain database). Therefore, we could argue that at least 15-50% of 
the He-shell material ejected from CCSNe was exposed to H ingestion. 

We explored the impact of H ingestion on the nitrogen inventory in the Galaxy. Intermediate- 


mass stars can largely account for the nitrogen abundance in the solar system (e.g., Kobayashi 


eTalllmTI ). Spectroscopic observations, on the other hand, have shown that massive stars are the 
most important suppliers of primary nitrogen in the early Galaxy ( Israelian et al.|[2004 Spite et al. 
2005). Recent GCE models underestimate the abundance of ^^N in the Sun by more than a factor 


of four (Kobayashi et al. 2011). Models of fast rotating massive stars can explain the observed 
production of primary ^^N (Chiappini et al. 2006), but not the missing ^^N (Meynet et al. 2006). 


This contrasts with observations of low ^^N/^®N ratios at high redshift, which suggests a massive 


star origin for ^®N (^^N/^®N=130±20 at 

z=0.89, 

Muller et al.|2006 

, while novae are a controversial 

source of ^®N at much later timescales ( 

Romano & Matteucci 

2003 

)■ 


We tested this scenario using the GCE model presented in Hughes et al. (2008), constructed 


with the GCE code GEtool (Eenner & Gibson 2003). We assume that the yields of all CCSNe 


with initial mass M < 20Mq carry the N isotopic ratio of model 25T-H (^'^N/^^N=20), and apply 
the prediction of model 25d for stars with initial mass M > 20Mq (^^N/^®N=434). Since N is a 
primary product of H ingestion and is not affected by the initial metallicity of the star, we used 
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the same solar yields for CCSNe of all metallicities. The GCE of the N isotopic ratio is shown 
in Fig.|^ in comparison with the baseline case. With the extreme assumptions mentioned above, 
we obtain an N isotopic ratio more than four times lower than in the Sun. Furthermore, our large 
production of and in the He shell of CCSNe changes the CCF results at low metallicities. 

We cannot rigorously quantify the impact of stellar yield uncertainties and CCF uncertainties 
in Fig.j^due, in part, to the limitations of current stellar models. Therefore, these results are meant 
to be illustrative rather than definitive. Hydrodynamic simulations of H-ingestion in massive stars 
are needed to produce more robust stellar yields to quantify the role that these events play in 
providing a source of nitrogen, and in particular of throughout the Calaxy. Major effects can 
also be expected for other isotopes, above all ^^C. 


5. Conclusions and final remarks 


We presented new explosive nucleosynthesis calculations for the He shell in CCSNe of massive 
stars. The main goal of this work is to study the nucleosythesis impact of the possible presence of 
H in the He shell, due to H ingestion before the SN shock reaches these layers. We explore a large 
range of explosion conditions (a factor of 3.3 change in the peak temperature of the SN shock, and 
a factor of 100 in the peak density), and of the H abundance (0.002%<Xh<1.2%) with simplistic 
ID SN simulations. 


Our calculations have identified the CCSN conditions to explain puzzling signatures in C- 
rich presolar grains. Although comparisons of these conditions to the observations can be used 
to constrain the nature of the supernova explosion, quantitative results require more detailed cal¬ 
culations. For example, the multidimensional structure generated by the H-ingestion cannot be 
properly represented by ID-hydrostatic models. The yield results are sensitive to the peak density 
and temperature and 3D-hydrodynamics simulations of the H ingestion are needed to quantitatively 
predict this pre-shock structure. With such detailed models, the features in the C-rich grains can 
be used to constrain the nature of the CCSN explosion and of SN-shock propagation producing 
these yields (e.g. Wongwathanarat et al. 2014). 


We explored the impact of our models on the CCF of the N isotopic ratio, showing that H 
ingestion in massive stars can be an important source of ^“^N and in particular of ^®N, affecting the 
solar N isotopic ratio. To investigate this idea in more detail, massive star simulations with mixing 
assumptions that allow H-ingestion events are needed. Such ID simulations need to be informed by 
3D-hydrodynamic simulations, and followed-up with more realistic CCSN simulations. For realistic 
CCF simulations, the yields of 10-12 models of high mass stars at 4 different metallicities (at least) 
are required. Our results indicate that a similar set of stellar yields from stellar models including 
H-ingestion events at different metallicities is needed. 
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M/Mq 



M/Mq 


Fig. 1. — Upper panel: Isotopic abundances in the He-shell ejecta of the 25Mq SN models 25T-H 
and 25T-H50 (thick and thin lines respectively). Shown are profiles for ^He, 
and the Si isotopes ^®’^^Si. The 0-rich 0/nova zone that formed in the 25T-H model is highlighted 
in grey (see text). Markers are used to identify different lines. Lower panel: we show the short-lived 
isotopes ^^Na, ^®A1, ^^Ti and ®°Fe in the CCSN He-shell ejecta. 
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Fig. 2.— Isotopic ratios and across the He/C zone predicted by the sets of 

models 25T (upper panel) and 25d (lower panel) are compared with measured data for single SiC 
grains (type AB, mainstream, nova, X and C) and LD graphites. Errors of grain data are shown 
(1-sigma). The solid vertical and horizontal lines indicate the solar isotopic ratios for C and N. 
For N we used Marty et al. (2011). The terrestrial N ratio (^'^N/^^N=272) is also shown (dashed 
horizontal line). 
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Fig. 3.— Upper Left: Isotopic ratios and from the 25T models are compared 

with those of C-rich grains (see Fig.j^. Continuous lines represents the C-rich He/C and C/Si zones. 
The dashed-line part of the curves in between indicates the 0/nova zone. Upper Right: Isotopic 
ratios and ^^C/^^C from the 25T models are compared with LD graphites. The 15Mq 

model by Pignatari et al. (2013a) is also reported (cyan left-pointing triangles). Lower Panels: 
Isotopic ratios ^®A1/^'^A1 and (Left Panel) and ^®Si/^®Si and ^^C/^^C (Right Panel) from 

the He/C zone of the 25d models are compared with AB and nova SiC grains. The Si isotopic 
ratios are in <5 notation (5(ratio)=(stellar ratio/solar ratio - l)xl000). For color/symbols of models 
and grains see Fig.[^ 
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Table 1: List of CCSN models. The He-shell SN ejecta are C-rich, excepting for the 0/nova zone 


the 25T model set. Details are given 

in the text. 



Xh in the He-shell 

before the SN shock 

1 .2% 

0.24% 

0 .12% 0.06% 

0.024% 

0.0024% 

CCSN models set d 
He-shell ejecta (M©) 

25d (P13) 

25d-H5 

25d-H10 25d-H20 

6.82-9.23 

25d-H50 

25d-H500 

CCSN models set T 
He-shell ejecta (M©) 

25T-H 

25T-H5 

25T-H10 25T-H20 

6.81-9.23 

25T-H50 

25T-H500 

0-rich 0/nova (M©) 

6.83-7.04 

6.84-7.0 

6.86-7.0 6.87-7.0 

- 

- 



Fig. 4.— GCE simulations of the N isotopic ratio using the CCSNe yields by |Woosley &: Weaver 


(|1995 ) (dashed line), and using their solar N yields modified according to the models 25d and 25T-H 


for all metallicities (continuous line). For intermediate-mass stars we used the yields by Karakas 


et al. (2010) 
























